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Real-time classical SU(iV) gauge theories at non-zero temperature contain linear divergences. We 
introduce counterterms for these divergences in the equations of motion in the continuum and on 
the lattice. These counterterms can be given in terms of auxiliary fields that satisfy local equations 
of motion. We present a lattice model with 6+ ID auxiliary fields that for IR-sensitive quantities 

^^ ' yields cut-off independent results to leading order in the coupling. Also an approximation with 

^\ , 5+1D auxiliary fields is discussed. 
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In recent years there has been a considerable interest in the calculation of the Chern-Simons diffusion rate in the 
symmetric phase of the electroweak theory. By the chiral anomaly this diffusion rate is related to the rate of baryon 
number non-conservation M and therefore it is an important parameter in electroweak baryogenesis scenarios; see 
for reviews e.g. Bpj . It was originally suggested by Grigoriev and Rubakov that such infrared sensitive quantities 
could be determined from a classical theory. Their idea was that field configurations at momenta k «T are expected 
q^ | to behave classically: 
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I. INTRODUCTION 
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0^ ' Since the diffusion of Chern-Simons number is dominated by soft fields with typical momenta k ~ g 2 T ', it was expected 

that it could be calculated in a classical theory [||]6). Later, it was argued by Arnold, Son and Yaffe that this 

c\ • picture is not complete, because soft field configurations are sensitive to (particle-like) thermal fluctuations with hard 

O ! momenta k ~ T, that are not correctly reproduced by the classical theory. These hard modes introduce Landau 

JL , damping and slow down the dynamics of nearly static magnetic fields at momenta k ~ g 2 T. Including these hard 

modes leads to an estimate of a typical time-scale t ~ (<7 4 T) _1 of these fields. Bodckcr S showed that also semi- hard 

modes at momentum scale gT contribute to leading order, and a even logarithmic correction t ~ (g 4 Tlog -) _1 to the 

time-scale is obtained. This logarithmic correction has been calculated in various approaches |9| 15 . 

X: 



17fl and if one takes the (continuum) cut-off A ~ T, the 
oops and with a proper matching this incorporates Landau 



To calculate the Chern-Simons diffusion rate in an (effective) classical theory with cut-off A, the corrections of the 
hard modes and semi-hard modes have to be incorporated. For the semi-hard modes this requires that the cut-off is 
large enough, A >> gT. To include the hard modes basically two approaches can be taken. One approach is that 
one uses the classical hard modes to mimic the effect of the hard modes in the quantum theory []16|] . These classical 
hard modes introduce linear divergences in the theory p 
classical hard modes generate the (classical) hard thermal 
damping in the correct way. 

A different approach is to include the hard thermal loops in the classical theory 18-22 . This can be done in an 
economic way by including an induced source in the equations of motion for the fields and to solve for the dynamics of 
the induced source in the presence of the classical fields. In that case one also has the problem of the linear divergences 
introduced by the classical hard modes. In the continuum these can be accounted for by adjusting the strength of the 
induced source ]2C| , |2^ |. 

On the lattice both these approaches have the problem that the linear divergences are not rotationally invariant 
p7[ . Therefore, the classical hard modes cannot correctly mimic the continuum hard thermal loop effects, like Landau 
damping, nor can these be to accounted for by adjusting the strength of the (continuum) induced source. The lattice 
spacing dependence of the classical theory has been studied in p4p5|] 
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The aim of this paper is to investigate the inclusion of counterterms for the linear divergences on a lattice. Via a 
continuum model (section |ll|) and a lattice model (section 111) suitable for perturbative calculations only, we arrive 



at two model s tha t may be used for non-perturbative calculations on a lattice. 



In section IV A the first non-perturbative model is discussed, it contains the exact counterterms for the linear 
divergences on a lattice with lattice spacing a^. The counterterms are generated by 6+ ID auxiliary fields that satisfy 
local equations of motion. The finite finite renormalization is chosen such that the theory is matched to a quantum 
lattice theory, with a smaller lattice spacing as- For very small couplings <?, the lattice spacing as can be chosen to 
be small compared to the inverse temperature and continuum results can be obtained. 



The second model is presented in section [VB it contains the counterterms in the approximation that the external 
momenta of a divergent diagram are strongly space- like (po << \p\). This model requires 5+ ID auxiliary fields and 
g is unrestricted. 

II. CONTINUUM 

In this section we review the HTL equations of motion in the continuum and show how the counterterms for the 
linear classical divergences may be obtained. It is convenient to use the local formulation of the HTL equation of 
Blaizot and Iancu [E6[. They showed that the HTL equation is equivalent to a set of linearized kinetic equations. The 
first equation describes classical gauge fields in the presence of an induced source 

[1^, ^(s)] =&*(*)> (2) 

with the covariant derivative D^ = d^ + igA^ and field strength F^ v = d^A v - d v ' A^ + ig[A», A v ], The induced 
current can be written as an integral over four-momentum K of the source density j v (K, x) 

jL & (x)=gJ^f{x,K). (3) 

The source density satisfies the equation 

[K M D»,f(x,K)] = 2gNK v KPF p(T (x)d^A(K) , (4) 

with N the number of colors and A(K) the free symmetric two-point correlation function of the fields that have been 
integrated out to obtain equations &|J). The factor 2 accounts for the spin degrees of freedom. Equations (0-0) are 
equivalent to the HTL's calculated with propagator A(K). 

The kinetic equations (0W) are derived under the assumptions that the theory is weakly coupled and that a; is a 
slowly varying variable, d << K ~ T '. Since we are interested in the divergences of the theory the loop momenta can 
be taken of the order of the cut-off k ~ A which is much larger than the fixed external momentum Q ~ d. 

In the continuum the two-point function takes the form 

A(K) = [Q(K ) - e(-K )]S(K 2 )N(K ) , (5) 

with N(Kq) the equilibrium distribution function of the fluctuations. The ^-function in (g|) forces the fluctuations 
on-shell with dispersion relation K° = ±fc, k = \k\. Inserting the propagator (pt) in (B), the solution for the source 
density can be written as 

j"(x, K) = 2NK U 5{K 2 ) \@(K )6N(x, k) - 0(-K o )5N(x, -ife)] . (6) 

The function SN(x, k) may be interpreted as a particle distribution function, because it satisfies the kinetic equation 

v"D M , SN(x, k)] = gv p F p0 N'(k) , (7) 

with the velocity v^ = (l,v) and v = k/k and N'(k) the derivative of the equilibrium distribution function to the 
energy. 

The induced source is proportional to the deviation from equilibrium of the density of particles in the plasma 



d 3 k 

(2^)3 



iind(z) = .9 / j^v'SNixj) , (8) 



with g the gauge coupling. 

Our task is to specify what we mean by the particles of the plasma. If we would consider QED, the particles would 
be the fermions and the equilibrium particle distribution function would be the Fermi-Dirac distribution function. 
In the SU(iV) gauge theory that we consider here, the particles are gauge field excitations. If the gauge fields in (g) 
and (fa) would be interpreted as mean fields then the particles can be identified with the thermal fluctuations and 
N(k) would be the Bose-Einstein distribution function riBE(k) |Pq| . However we consider the equations (§-0) as the 
equations of motion of a classical statistical theory, where still a thermal average over the initial fields has to be taken. 
The particles in this case are the quantum fluctuations with distribution function 

N(k) = n BE (k) - T/k , (9) 

as can be derived from separating off the classical fields from the path integration and calculating the two-point 
function of the remaining fields |27|] which appears in (Q) . 

The equations (&R) may be simplified by introducing the field 

W(x, v) = -5N(x, k) {gN'ik))" 1 . (10) 

It follows from (H) that W is independent of the energy k, assuming this is consistent with the initial data. The 
kinetic equation (|7|) in terms of this field reads 

[v»D^W(x,v)]=v-E(x) . (11) 

The induced current (H) can be written as 



-2 9 'NJ^N'(k)J^W(x,v), 



jy nd (x) = -2g*Nj ^k 2 N'(k)J —v v W{x,$) , (12) 

where the angular integration is over the direction of v. As usual the fc-integration decouples from the angular 
integration in the HTL approximation. 

Due to the subtraction in N(k) (|9|), the integration over k introduces a linear divergence. This divergence acts as a 
counterterm for the linear divergences in the classical theory p8j. Let us argue why this is correct. The argument used 
in J28| is one of consistency: since the full quantum correlation functions cannot contain classical linear divergences, 
the quantum corrections to the classical equations should provide the correct counterterms. Another argument is 
that integrating out the classical fluctuations around a given background field generates the classical induced source 
in the effective equations for the background field. Since the classical induced source (in the HTL approximation) is 
precisely what is subtracted in (|l2]), the equations of motion for the background field do not contain linear divergent 
terms. The logarithmic divergences of the classical theory are more complicated and are not considered here. 

If we use dimensional regularization, the subtraction does not contribute, since linear divergences are set equal to 
zero. We then have 



Jind ( x ) 



Sc^I^Wm, (13) 



with the plasmon frequency uA = g 2 NT 2 /9. 

As a different regularization scheme that makes the linear divergence explicit, we may restrict the subtraction to 
momenta k < A: 

N A (k) = n BE (k)-^Q(A-k). (14) 

The result of the /c-integration in the induced source yields then 

2^ 
where we see an explicit linear divergence. The inclusion of a linear dependence on the cut-off in the plasmon frequency 



2 9 2 N J ^k 2 N' A (k) = 3u& - V^VAT , (15) 



as in (|15|) was proposed in 20 23J . The derivation presented here and in more detail in |28|] shows how such a cut-off 
dependence arises in a continuum theory as a consequence of the subtraction in the distribution function (0) for the 
quantum fluctuations. 

From the point of view of renormalization the situation is quite remarkable: there is an infinite set of linearly 
divergent diagrams that are all non-local (or momentum-dependent) but the renormalized equations of motion can 
be brought in a local form and the linear divergence is introduced in just one parameter. Note also that the usual 
HTL's are just a finite renormalization. 



III. PERTURBATIVE RENORMALIZATION ON A LATTICE 

Before turning to the HTL equations of motion, we shortly review the static classical theory on a lattice, as far as 
the linear divergences are concerned. The only linear divergence in the static theory occurs in the Debye mass. A 
counterterm for this divergence may be introduced in the mass of the temporal gauge field |29| , |30| 

8m = m D — ra cl , (16) 

with the continuum HTL contribution 

m 2 D = -2g 2 N J 4^n'(k) = \g 2 NT 2 , (17) 

and the classical mass (for a simple cubic lattice with lattice spacing a) 

m 2 cl = -2g 2 N J -^- n '(tip) « OMg^NTa' 1 , (18) 



where p is restricted to the first Brillouin zone \pi\ < ir/a and the energy tip is 

[atf (M) W (M) W (5*2)] . (19) 



p a 2 

The mass (|l8|) is the linear divergent contribution to the Debye mass on the lattice, its subtraction in (|16|) ensures 
that no linear divergences are present in the static theory with the mass counterterm included. The continuum HTL 
contribution ( |17| ) to the counterterm mass ( J16J ) provides the finite renormalization, which is chosen such the a — > 
limit yields the continuum result. 

We want to extend this approach to the real-time classical theory. We consider again the the equation for the 
gauge fields (H), but with a lattice regularization; time is continuous but space is a simple cubic lattice with lattice 
spacing a. Similar to the mass counterterm ( |lq ) in the static case, we want the source to contain a continuum HTL 

contribution with a classical lattice contribution subtracted. To this end, we introduce two particle densities SN(x, k) 
and SN ct (x,p) for particles with energies E% = \k\ and tip respectively. The idea is that the particle density SN ct 
generates the counterterms for the linear divergences and SN generates the correct finite renormalization , the "good" 
HTL contributions. Then these particle densities satisfy the equations 

v"D„,6N(x,k)\ =gv p Fp a {x)n' BE {k), (20) 

[t&Dr, 5N ct (x,p)} = gv^F p0 (x)n' cl (tip) . (21) 

Here and in the following, we persist in using a continuum notation although we use a lattice regularization for the 
UV-divergences. The kinetic equation for SN(x, k) is eq.(R) with the Bose-Einstein distribution function as equilibrium 
distribution function. In the kinetic equation for SN ct (x,p) the velocity on the lattice is vf^ t = (l,wi a t) with |lq ] 

<t = d P, ^p = -FT sin(aft) , (22) 

til L iff 

and |wi a t| 7^ 1 in general. The induced current in the classical equations of motion for the gauge fields then read 

j? nd (x) = 2gN J ^L v »6N(x, k) - 2gN J ^v? at 6N ct (x,p) . (23) 

Here the integration over p is restricted to the first Brillouin zone \pi\ < n/a. 

As in the continuum, a field W(x, v) may be defined that satisfies (the lattice version of) equation ( |ll|). Since 
the lattice velocity (E3) is not restricted to the speed of light, we have to allow for general velocities v in (|ll|). The 
induced current (g3|) reads 



J fV^(*, V) - 2g 2 NTa-' J ^{lZ 2 v ^W( Xl t/ lat ) 



j? nd (x) = 3^ J —v v W(x, v) - 2g 2 NTa- 1 J ^VZ 2 v^ t W (x , v lat ) . (24) 

with the dimensionless quantities pi = api, tip = atip and the integration restricted to \pj\ < tt. The first term on the 
r.h.s. of (p4|) is the continuum contribution for which the ^-integration decouples, as in (|l2|), and has been performed. 



In the second term on the r.h.s. of ( |24| ) the integration cannot be simplified since the velocity not only depends on 
the direction of the momentum p, but also on its magnitude. For the calculation of the continuum contribution to 
the induced current a field W(x, v) depending on the direction of v only is sufficient, howev er the lattice contribution 



requires fields that depend also on the magnitude of the velocity \v\ a t\ < 1. In section IV B we will study the question 



if, for the calculation of the Chern-Simons diffusion rate, we may approximate the induced current with fields that 
only depend on the direction of the velocity. 

We note that the induced current ( p4[ ) is covariantly conserved \D^ , j? nd ] = 0. 

Just as the usual HTL equations, the equations (||), (JTTJ) and ( |24| ) (or equivalently (g), (|20|), ( [2l|) and (|23|)) define 
a perturbation theory. Taking retarded initial conditions the retarded propagator (and higher-order retarded vertex 
functions) can be obtained, as in |2(|. The classical KMS condition then fixes the entire propagator, including its 
thermal part pjfl . Using perturbation theory we may verify that also the time-dependent counterterms are correct, 
we calculate the retarded propagator to one-loop order. In a general gauge it takes the form 



D7(Q) = [g^Q 2 - Q^Q V + w + Ki(Q) + «C d (Q)] * 



(25) 



with F^ the gauge fixing vector and II^" the classical self-energy and ^n^ the counterterm self-energy introduced in 
the induced source (|^). The classical self-energy to one-loop order reads ]17| , p4| 



Kl(Q) = ^Na- 1 
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_JMO^0 



U iat V lat<?0 



q a + ie - v\ &t • q 



(26) 



at this order the classical self-energy contains no contribution from the induced source. The linearized induced source 



iind(^) 



^x'5-U.r A {x,x')A v {x') 



ind \ 



defines the retarded self energy 



with the continuum HTL self-energy 



(27) 



(28) 



(29) 



Inserting the counterterm self-energy ( p8|) in the propagator (J25|) , we see that the linear divergent classical self-energy 
and the linear divergent part of the counterterm self-energy cancel. The resulting self-energy in the propagator (Eq) 
is the correct (continuum) HTL self-energy. Finally we note that in the static limit the self-energy (|28|) reduces to 
the counterterm mass (fig), as it should. 

Unfortunately the system (p|) , (Jll]) , ( p4| ) is unsuitable for numerical implementation [Q. This can be seen from the 
conserved energy of the system 
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htl(Q) ~ 3^ pl 
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{E b Y + (B b y + 3lu 



pi 



— W b (x,v)W b (x,v) - 2g 2 NTa~ 1 
4tt 



d3 P A-2tt... - .,,/-, 
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T VLZ 2 W\x,v ut )W b {x : v ut ) 



(30) 



with B the chromo-magnetic field and b the adjoint index. The energy is unbounded from below and this means that 
the system is unstable. In an ideal world the effect of the counterterm particle density is precisely compensated by 
the hard modes of the classical gauge fields, however in practice the evolution of the particle density and the hard 
modes will differ, which means that after some time the (wrong) effect of the counterterm particle density is no longer 
compensated by the hard modes, and the fields will (exponentially) blow up. 



IV. TWO MODELS WITH A BOUNDED ENERGY 



A. Model with lattice dispersion relation 



The goal is to obtain a model that is defined on the lattice, that is stable and that can be used to calculate 
IR-sensitive real-time properties of a non-abelian plasma without linear divergences. Such a model should meet the 



following three requirements: 

1) In the small lattice spacing limit the continuum HTL equations of motion should be obtained. 

2) Counterterms for the linear divergences (on the lattice) should be included. 

3) The energy must be bounded from below. 

The model considered in the previous section failed to have bounded energy. To obtain a model with a bounded 
energy one can consider a model where the modes inducing the finite renormalization have the same dispersion 
relation as the counterterm modes. In this section we focus on a model where both the counterterm modes and the 
modes generating the finite renormalization satisfy a lattice dispersion relation. The other possibility of enforcing the 
continuum dispersion relation on the counterterm modes is considered in the next section. 

To obtain HTL equations where the both types of modes satisfy a lattice dispersion relation, we do not match to a 
continuum quantum theory as in the previous section, but to a quantum theory on the lattice, with a (small) lattice 
spacing as- The HTL equations (in the Aq — gauge) may then be written as 



[D fl ,F^(x)]=Jl' nd (x) = 2gN 



d 3 p „ 



d t SN(x,p) ~ [«iat • D,5N(x,p)} = -gvut ■ E{x) d {l N(n p ) , 



v^SN(x,p) 



(31) 
(32) 



with x = (t,x), where the time t is continuous and the position x is an element of a cubic lattice with (large) lattice 
spacing a^. The dimensionless momentum p is restricted to the first Brillouin zone \f>i\ < w, the dimcnsionless energy 
2\7Si sm (P«/2) 2 and the velocity vl 



i Lp 



"lat 



8- Q- 



The lattice spacing has been scaled out of the above equations and enters only in the equilibrium distribution 
function N. The distribution function N should contain a contribution that, after solving (|32j), generates the quantum 
HTL source and a contribution that generates the counterterms for the classical divergences. The important step is 
now to allow for different lattice spacings a^, as in the the different parts of the equilibrium distribution function 



N{Qp) 



-2„S 



n B E(^p)~a L n cl (fl p ) , 



(33) 



with 



n BEV^f>) 

n^fip) 



1 



e flp/(a s T) _ I ' 

Ta L 



- ' n 



(34) 



with T the temperature of the system 
To see that the model @ and (P 
the field 



contains the counterterms for the linear divergences it is useful to introduce 



W{x,p) = 5N{x,p)/(-gN\h p ) 



(35) 



which satisfies the equation 



d t W(x,p)- wi a t -D,W(x,p) =vut-E(x) 



(36) 



The source can be split into a part generating the finite quantum HTL source and a part subtracting the linear 
divergent classical source 



7ind Jfin 7ct 



(37) 



In terms of the field W these sources read 



;l = 2^/^.< •■' 



f ct = 2g*N I d3pL ^ ' 



7^T3^ llt n' BE {ns)W(x 1 psas) 
-V\ a ,tn'dp. L )W{x,pLa L ) 



(2tt) 3 



(38) 
(39) 



with ps = a s p, Qs = a s ^p an< ^ similar for p^, fl^. Both sources (pq) and (p9J) are covariantly conserved 



Written in dimensionfull quantities we recognize the source j c t ( J39J ) as the classical HTL source on a lattice with 
lattice spacing ah- The difference with the perturbative model of the previous section is the choice of the finite 
renormalization. The source j nn (pq) is the quantum HTL source on a lattice with lattice spacing as- To extract 
continuum results from this model we should require a^ >> T . Also a^ cannot be too large, since the relevant field 
configurations for the sphaleron rate have size (<7 2 T) _1 , we should at least require a^ >> g 2 T. However as Bddeker 
H has shown modes of spatial size (gT) -1 give corrections of 0(1), to take these corrections into account requires a 
smaller lattice spacing a^ >> gT 

To ensure the stability of the model ([u]) and (|3^) we demand that the energy, 



E = 
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{E b Y + (B b ) 



2N 
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J2nf 



6N b (x,p)6N b (x,p)/N'(Cl p ) 



is bounded from below. This leads to the requirement 



N'(np) > 



(40) 



(41) 



For p = 0, this requirement implies 05 < a L , which is in accordance with the general idea that the classical theory is 
matched to a quantum theory with a smaller lattice spacing. 

The function — N'(Qp), with as < ol, decreases from plus infinity at flp = 0, to its minimum below zero, after 
which it increases and asymptotically reaches zero. The maximum value of the dimensionless energy is flp = 2y3- 
Demanding that 



- AT'(2\/3) > 



(42) 



together with as < a L is sufficient for (|4l]) to hold for any p. In this way, we obtain a maximum value for a^ 1 given 
the ratio a^j as- In table H some results are listed. We see that in order to obtain continuum-like HTL contributions, 
the ratio a L jas should be very (exponentially) large. 

Since we want a^ >> gT, the coupling coupling g should be chosen extremely small. For instance, if we fix 
a], 1 = 2.59T, then stability requires a L /a s > 100, so Oj} < 2.59 10~ 2 T and g « 2.59 10~ 2 . 

TABLE I. The maximum value of a^ 1 /T given the ration aL/as from the requirement that the energy is bounded from 
below. 



0-L/as 


1.1 


1.5 


2 


5 


10 


20 


25 


50 


100 


1000 


max(a^ i )/T 


0.31 


0.64 


0.86 


1.36 


1.68 


1.97 


2.06 


2.33 


2.59 


3.42 



To complete the model we specify the average over the initial fields, that has to be taken to calculate classical 
correlation functions. The initial conditions are 



E(t in ,x) = E- m {x) , 

A(t\ ni x) = Ai n [x) , 

8N(t- lni x 7 p) = SN in (x,p) . 

Following pf3J33]] , we include the auxiliary field in the average over initial fields 



Z = J VE in T)A- m mN m 6 (G in ) cxp (-/3E) , 



(43) 



(44) 



with Gi n = Gauss' law at the initial time. We have verified that the phase space measure is invariant under time 
evolution. 

We may conclude that the model given by ([32|) , (J43|) , (|44|) describes a real-time quantum theory on lattice to leading 
order in g without cut-off dependence. For very small coupling the lattice spacing as of the quantum theory can be 
taken small compared to the inverse temperature and the model describes a quantum continuum theory. It can be 
used to calculate Hi-sensitive quantities without lattice spacing dependence. Consider for instance the Chern-Simons 
diffusion rate, which in the small coupling limit takes the form pi 



r = 

A calculation of T with the model given by ( |32| ) 
for both coefficients n\ and k 2 , for a^ 1 >> gT. 



Kx log - + K 2 
9 



10rp4 



g w r 



(45) 



(43) 



at fixed lattice spacing as, gives an a^-independent result 



B. Model with a continuum dispersion relation 

The other approach that we want to investigate is a model where we enforce the continuum dispersion relation on 
the counterterm modes. Such a model has the advantage that instead of a 6+1D field SN a 5+1D auxiliary field 
W(x,v\ a t), that depends only on the direction of the velocity V\ a t — Qat/I^latl) can be used. The counterterms that 
we obtain in this model are not exact, but for the calculation of the Chern-Simons diffusion rate the model provides 
a reasonable approximation. 

The model that we consider is given by the replacement of the induced source (O) by the expression 

f app {x) = 3c£ J ^VW{x, v) - 2g 2 NTa- 1 J -^QZ 2 \v lat K t W(x, v iat ) , (46) 

with v" at = (1, fiiat)< We choose this expression since it reproduces the induced vector current for a field configuration 
with W(x, -Oiat) = W(x, wiat), and the vector current is essential in the dynamics of the soft fields. The density is then 
determined by requiring current conservation [Z) M , j app ] =0. As a consequence the induced density j® nd in ( gj ) is not 
correctly reproduced by the density j app . This can be easily understood, changing the velocity of the particles and 
requiring current conservation either the vector current or the density can remain unaltered. The expression (J46J) is 
the lattice equivalent of the approximation for the induced source in EG] . 
We may also write (flq ) as 

JappW = / ^m 2 (v)v u W(x, v) , (47) 

with the velocity dependent mass 

m 2 (v) = 3o& - 2g 2 NTa- 1 J ^n-J\v lat \6 s (v~v lat ) . (48) 

The second term contains a linear divergence in the direction v — (1, 1, l)/v3 p2] and logarithmic divergences in 
directions v = (1, 1, s)/y2 + s 2 with — 1 < s < 1 (and directions related by symmetry). Therefore the mass and the 
energy are not strictly positive. To obtain a bounded energy some averaging over the direction of the velocity v should 
be performed. This can be achieved by expanding the field W(x,v) in spherical harmonics 

W(x,v) =J2 Wlm W Ylm< $) , (49) 

I in 

and keeping a finite number terms. The induced source can then be written as 

3^d(x) = ^2ar m W, m (x), (50) 



liri 



with coefficients 



aL = /^m 2 («KY im (iT). (51) 

Given the lattice spacing a, the requirement that the energy is bounded from below puts an upper bound Z max on 
allowed values of I, It was found in J23] that the Chern-Simons diffusion rate is insensitive to l m3X for even l max . In 
the following we will therefore focus on the approximation made in ([46[). 

As was already mentioned the approximation (Eq) changes the charge density. For instance for the coefficient a§ , 
we have 

a° =m% - 2g 2 NTa- 1 J ^Lfl~J\v ut \ . (52) 



Comparing ( |52| ) with (18), we see that the expression ( |46[ ) does not correctly reproduce the counterterm for the Debye 
mass. This implies that the current is not suitable to describe the behaviour of fields at length scale (gT)^ 1 . 

To see whether the approximation (Efl) is valid for fields at the length scale (g 2 T)~ 1 , we consider the spatial 
components of the counterterm self-energy generated by the source (ffq) (for Z max — > oo) 



nSGto, q) = n^ TL ( go , s) - ng p ( qo , <f) , (53) 

with 

H% p ( qo ,q) = 2g 2 NTa^ f^-^\vu,\ §^f . , (54) 



which should be compared with the classical self-energy (|26[ ). 

It is important to realize that the relevant fields for the Chern-Simons diffusion rate we are interested in, have typical 
momenta of order qo ~ g 4 T, q ~ g 2 T. For the gauge fields relevant for the Chern-Simons diffusion rate qo << \<j\ and 
neglecting q in the denominator of the counterterm (p4h and the classical self-energy (Eq), we note that these are 
equal and that they cancel. For these fields the effective theory is finite and reproduces the HTL contributions. 

However it was realized by Bodeker that interactions between semi-hard and soft fields give corrections to the dy- 
namics of the soft fields that are not suppressed by powers of g. On the contrary, even log(l/<?) enhanced contributions 
arise. The counterterms in the approximated source (Era) and the classical HTL's do not cancel for the semi-hard 
modes (with momenta qo, q ~ gT), therefore the semi- hard modes are sensitive to the cut-off a -1 . 

The leading log contribution arises from the IR-sensitive part of the contribution of the semi-hard modes with 
momenta fco << k ~ /i and /j, ~ g 2 T an IR cut-off. For these momenta the approximation is correct to leading order, 
therefore a calculation of the Chern-Simons diffusion rate with approximation (|4q ) produces the correct leading log 
contribution, the coefficient K\ in (fiq ) is independent of the lattice spacing. 

The 0(1) correction from the semi-hard modes does depend on the cut-off. An estimate of the cut-off dependence 
can be obtained from a comparison of the classical HTL self-energy (E6I) with the counterterm (p4j). To be explicit, 
we compare the diagonal components at zero spatial momentum 

n*t(go,<f= 0) = 2g 2 NTa- 1 j ' -0-{l~ 2 \v lat \ 2 = 0.26 5 2 ATa" 1 , (55) 

K pp (q ,q = 0) = 2g 2 NTa- 1 J -0-ClZ 2 \v lat \ = O^NTa' 1 . (56) 

Comparing the difference between (45) and (46) with the HTL self-energy at zero spatial momentum TTfj TL (q Qj q = 
0) = 3ojp] = g 2 T 2 /3, we obtain an estimate for the maximal error of about 25% for a -1 = T/h. However, the 
semi-hard modes that give the 0(1) correction have space-like momenta qo < \q\ Ej]. For these modes we expect (p3J) 
to be a better approximation of the classical self-energy (£6|) . 

Besides the mismatch between classical HTL's and the counterterms from (^), the lattice spacing dependence of 
«2 depends on the magnitude of the 0(1) correction from the semi-hard modes. Especially when the soft modes 
dominate the contribution to k,i this model is suitable for a calculation of the Chern-Simons diffusion rate. 

V. CONCLUSION 

In this paper, we studied the linear divergences in classical SU(iV) gauge theories at finite temperature. Countert- 
erms for these divergences can be incorporated in an (induced) source. Although the divergences are non-local the 
equations of motion including these counterterms can be given in a local form by introducing auxiliary fields. In the 
continuum these auxiliary fields depend only on the direction of the velocity v, whereas on the lattice they necessarily 
also depend on the magnitude of the velocity. We have presented two lattice models with a bounded energy. The 
first describes a real-time quantum lattice theory with a small lattice spacing as- The requirement that the energy is 
bounded presents a lower bound on as given the lattice spacing ai, of the classical model. It is shown that for very 
small coupling g the lattice spacing as can be taken small compared to the inverse temperature, which implies that 
continuum results can be obtained. 

In the second model we argued that the restriction to auxiliary fields depending on the direction of the velocity 
allows for a reasonable approximation (|46|) for the calculation of quantities dominated by fields with momenta (qo, q) ~ 
(g 4 T,g 2 T), such as the Chern-Simons diffusion rate. 
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